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ABSTRACT 

Surface shedding of an asteroid is a failure mode where surface materials fly 
off due to strong centrifugal forces beyond the critical spin period, while the 
internal structure does not deform significantly. This paper proposes a possible 
structure of an asteroid interior that leads to such surface shedding due to rapid 
rotation rates. A rubble pile asteroid is modeled as a spheroid composed of a 
surface shell and a concentric internal core, the entire assembly called the test 
body. The test body is assumed to be uniformly rotating around a constant 
rotation axis. We also assume that while the bulk density and the friction angle 
are constant, the cohesion of the surface shell is different from that of the internal 
core. First, developing an analytical model based on limit analysis, we provide 
the upper and lower bounds for the actual surface shedding condition. Second, 
we use a Soft-Sphere Discrete Element Method (SSDEM) to study dynamical 
deformation of the test body due to a quasi-static spin-up. In this paper we 
show the consistency of both approaches. Additionally, the SSDEM simulations 
show that the initial failure always occurs locally and not globally. In addition, 
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as the core becomes larger, the size of lofted components becomes smaller. These 
results imply that if there is a strong enough core in a progenitor body, surface 
shedding is the most likely failure mode. 

Subject headings: minor planets, asteroids: general 


Introduction 


The Yarkovsky-O’Keefe-Radzievskii-Paddack effect, so-called the YORP effect, play a 
crucial role in prov idin g asymm etric asteroids with rotational torque due to the solar radi¬ 
ation pressure ( Rub i ncam 120001) . This implies that asteroids wit h sever al hundred meters 
in diameter may be able to gradually spin up to their spin limits (iRubincaml 200011 . Obser¬ 
vational surveys have shown that there is a spin period threshold of ~ 2.3 hr for asteroids 
larger than a few hundred meters in size (iPravec et al.l 20071). Recent observations de¬ 


tected spin-up/down of asteroids due to th is effect flLowrv et al.l 120071: iTavlor et al.l 120071: 
Kaasalainen et ah 2007; Durech et ah 2008]). These reports imply that spin-up processes 
could cause asteroid structures and shapes to evolve over their life time. 

Some recently observed disruption events of active asteroids are considered to have re¬ 


sulted from ro tational instability; how ever there are o 


events as well. IJewitt et al.l (120100 and ISnodgrass et ah 


h er pr oposed mechanisms for these 
2010f) both independently reported 


the disruption event of asteroid P/2010 A2. IJewitt et al.l (120101) concluded th at it was l ikely 
that t he formation of the debris tail resulted from rotational instability, while ISnodgrass et ah 


( 2Qlol) noted that a collisional impact is a likely event driving such events. Jewitt et ah 


( 20141 ) described the recent breakup event o f P/2013 R3, suggesting that this event came 
from rotational instability. Hsieh et ah (200 4 . 2010) reported that the dust tail of asteroid 
133P/Elst-Pizarro, spinning with a spin period of 3.471 hr, may come from seasonal activ¬ 
ity and rotational instability. Sheppard & Trujillo (120151 ) also observed a thin dust tail of 
asteroid (62412) 2000 SY178 whose spin period is 3.33 hr from March 28 to May 2 in 2014. 
Although the activities of these bodies could be driven by sublimation of water ice that 
might be buried in the bodies, the slow particles in their dust tails also imply th at rota tional 
instability may also play a role in mass ejection from the surface. IJewitt et ah (120131120151 ) 
conducted observations of P/2013 P5 (PANSTARRS) with the Hubble Space Telescope to 
investigate its episodic du s t ejec tion, implying that this event mig ht result from rotational 
instability. lHainaut et al.l (120141) and lHainaut fe Snodgrassl (120141) proposed an alternative 
scenario of this event in which the dust tail might be generated by the soft contact in a 
binary system gently rubbing each other. 


Over the last decade, there have been many important analytical (Ho lsap ple 200ll l2004 
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20071. l2010h and numerical flWalsh et al.ll2008l. 120121: ISanchez fc Scheeresll2012l) investigations 
of failure conditions of asteroids. Of specific relevance for this study, it has been found that 
for a homogeneous strength distribution in an astero i d, t he internal structure is more sensitive 
to failure than the surface region ( IHirabavashi fe Scheeresl 20151) . Another active research 


area is an understanding of mass ejection from a fast rotating body due to centrifugal forces. 
A popular explanation for the origin of mass ejection is that surfa ce materials flow down to 
the equator, and if the rotation is fast enough, they fly away. Guibout & Sch eer esl (120011 ) 
investigated the surface stability due to the shape and spin period of a uniformly rotating, 
self-gravitating ellipsoid and found that based on the body configuration, the surface ma¬ 
terials may acc umulate at d iff erent ends. Using a Hard-Sphere Discrete Element Method 
(HSDEM), W alsh et al. ( 20081 . 1201 21) showed the formation of a binary system due to m ass 
ejectio n fro m a cohesionless sp heroid. By considering a constant angle of repose, iMinton 
( 2008 ) and Harris et ah (12009 ) numerically obtained equilibrium shapes of a uniformly ro¬ 


tating symmetric body and obtained a body shape similar to 1999 KW4 Alpha. Schee res 


( 20151 ) analytically investigated granular surface flows on rotating asteroids and their equi¬ 
librium shapes, demonstrating the formati on of unique surface slopes of oblate shapes such 
as 1999 KW4 Alpha and 2008 EV5. Also, IHirabavash i fe Scheeres (2 0141 ) compared surface 
mass ejection and structural failure due to rapid rotation to determine a possible failure 
mode of (216) Kleopatra. 

However, the mechanism behind surfaces mass ejection is still poorly understood, and 
thus we will try to s hed s ome light on it in this paper. Developing a two-bulk-density-layer 
model, Hirabavashi (|2014 ) showed that a hi gh-density core makes the whole body structurally 
strong. Assuming a homogenous structure, Hirabavashi fe Scheeresl ( 2015 ) conducted plastic 
finite element analysis to explain that the failure mode of asteroid (29075) 1950 DA is plastic 
deformation of the internal region. These studies imply that in order for the surface region 
of an asteroid to fail Erst at an elevated spin rate, its interior must be stronger than a 
surface layer. Based on these studies, we propose one way to create surface shedding from 
a rubble pile asteroid. Here, we introduce a simplified model to discuss it. Our rubble pile 
asteroid is modeled as a self-gravitating spheroid composed of two cohesive strength layers, 
a spherical shell and an internal core. The bulk density and the friction angles are assumed 
to be constant to simplify our discussion. We note that cohesion may depend on porosity 
(Lambe & Whitman 19691 ). but we will leave consideration of cohesion and porosity for a 
future study. 

In the following discussion, we will use the term “sur face shedd ing.” So, before we do, 


we will try to give it a proper definition. Hirabavashi fe Scheeresl (120141 ) introduced this 
term to explain the mode that particles just resting on the surface are shed due to strong 
centrifugal forces. However, this definition is still ambiguous. First, if particles on the sur¬ 
face interact mechanically with other particles by cohesion, it is difficult to determine the 
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lofting condition by simple balance between gravitational and centrifugal forces. Second, 
if cohesion is uniformly distributed over the whole volume, the internal structure deforms 
before the force balance point, i.e., the dynam ical equilibrium point, reaches the surface 
fjHolsapplel l2010l; iHirabavashi fe Scheeresl 120151) . To define this mode clearly, we describe 
“surface shedding” as a mode in which while the internal structure is below the yield condi¬ 
tion, the surface region fails structurally and particles in this region are shed due to strong 
centrifugal forces. We hypothesize that this mode may p lay a role in the disruptio n events of 
133P/Elst-Pizarro (Hsieh et al. 2004 . 2010 b P/2013 P5 ( Jewitt et al. 2013 . 2015 ) and 62412 
(Sheppard fe Truiillo[2015h 


2. Analytical Modeling of Mechanical Failure 
2.1. Internal Core Model 

We suppose that the test body considered here is a uniformly rotating, self-gravitating 
spheroid that has a concentric, spherical core under the surface shell, later known as the 
internal core (Figured]). It is assumed that these layers have different cohesion, while their 
bulk density and friction angle are the same over the whole volume. Surface shedding is 
induced not by the failure of the internal core but by that of the surface shell. If the test 
body is spinning fast enough, surface particles are shed from the test body; otherwise, they 
just move over the surface. 

The bulk density, the total radius and the gravitational constant are denoted as p , R 
and G, respectively. We normalize lengths, body forces, spin rates and stress tensors by 
R, 7T pGR, y/irpG and np 2 GR 2 , respectively. With this normalization, the sphere radius is 
denoted as 1 and the radius of the internal core is defined as Rb, so the thickness of the 
surface shell is given as 1 — Rf,. The right plot in Figure Q] shows the spherical coordinate 
system (r, 0,</>), but we also use the cartesian coordinate system ( x,y,z ) as well. The spin 
axis is assumed to be constant along the z axis. We consider that at a given spin rate, if 
the surface shell already fails structurally, the cohesion between the particles there is broken 
and they can move freely. For this case, the particle is lofted from the surface due to the 
centrifugal force at a spin rate of a/4/3 ~ 1.15, while it energetically escapes to infinity at 
that of ^/8/3 ~ 1.63. 


2.2. Failure Conditions of the Two Layers 

In the analytical model developed here, the behavior of materials is assumed to be 
elastic-perfectly plastic and to follow an associated flow rule. The yield condition of a 
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material is assumed to be characterized by the Mohr-Coulomb yield criterion, a pressure- 
shear dependent yield criterion, which is given as 


(< 7 i — < 73 ) sec (j) + (< 7 i + < 73 ) tan (j) < 2 Y. 


( 1 ) 


where a t (i = 1, 2, 3) is the principal stress component, 0 is a friction angle and Y is cohesion. 
Since we consider a sphere, the stress state is always < 7 i = < 7 2 > 03 . 


The yield conditions of the surface shell and the internal core are calculated by using an 
elastic stress solution. Solving the elastic equations, including the equilibrium equations, the 
constitutive equations (Hooke’s law), the strain-displacement eq uations and proper boundary 
conditions (zero-traction condition), we obtain the solution as (iDobrovolsk is 19821 ) 


&i(l — r 2 ) — ki 3 x 2 — k u z 2 , 

(2) 

h{ 1 - r 2 ) - k 13 y 2 - k 15 z 2 , 

(3) 

k 9 (l - r 2 ) - k 14 x 2 - k 15 y 2 , 

(4) 

hs xy, 

(5) 

k^xz, 

(6) 

ki 5 yz, 

(7) 


where r 2 = x 2 + y 2 + z 2 and the coefficients, ks, are defined as follows: 

2b x (— 12 — v + 5iA) + b z ( 3 — 6u — 5iA) 

1 ~~ 10(-l + z/)(7 + 5z/) ’ 

26 x (3 — 6u — 5u 2 ) + b z (— 27 + Au + 15i^ 2 ) 

9 “ 10(-l + z/)(7 + 5z/) ’ 

2 b x (y —4 + 3^ + 5^ 2 ) + b z ( 1 + 3^) 

13 = 5(-1 + i/)(7 + 5i/) ’ 

_ &x(-3 + 61 / + 5i/ 2 ) + b z (-A + 3z/ + 5^ 2 ) 

14 ““ 5(-l + z/)(7 + 5z/) ' 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Also, k§ = k\ and ki A = ki 5 . v is Poisson’s ratio. The coefficients of the body force 
components along the x, y and z axes, denoted as b y and b z , respectively, are given as 


b x = b 


= UJ 
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( 12 ) 


To determine lower and upper bound condi tions for the actual failure, we apply the 
upper and lower bound theorems in limit analysis ( Chen fc Han 1 988 : Chen & Liu 19901). At 
a lower bound condition, the target volume should not fail structurally. At an upper bound 
condition, on the other hand, the volume should fail. Thus, the actual failure condition 
should always be between these conditions. 
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The lower bound theorem describes that when (1) an equilibrium stress solution is found 
and (2) there is a small element whose stress reaches its yield condition, the body may not 
widely fail. This means that the yield condition of a small element does not guarantee 

further propagation of failure regions. Therefore, we say that this condition is below the 

actual failure condition. Practically, this condition is given as the yield condition of the 
stress at the most sensitive element. Because the test body is spherical, such an element 
always appears on the equatorial plane. 

For the lower bound condition of the surface shell, the boundary between the surface 
shell and the internal core becomes the most sensitive to failure. Considering y = z = 0, we 
obtain the principal stress components at this location as 

cr e i = fci(l-r 2 ), (13) 

cr e2 = h( 1 - r 2 ) - k 13 r 2 , (14) 

cr e3 = kg( 1 - r 2 ) - k 13 r 2 . (15) 


It is not always satisfied that a e \ > cr e2 > cr e 3 . Thus, to have the correct order, we calculate 
a es i = max(<7 e i, cr e2 , a e3 ) and a es3 = mm(a e i,a e 2 ,a e3 ). On the other hand, for the lower 
bound condition of the internal core, since the center is always the most sensitive location, 
we obtain the principal stress components as 

(T e cl k |. < 7 ec2 k 1 , a ec3 kg. ( 16 ) 


The lower bound condition provides the highest approximated value of cohesion for the 
actual failure condition, later known as the highest cohesion. Based on Equation (0Q), this is 
expressed as 

1 1 

^efc ^ &ek(l) S6C (f) T “ (o"efcl T (Tekz) tan 0, (17) 


where k is either s or c. k = s indicates cohesion of the surface shell, while k = c provides 
that of the internal core. 


F or the u pper bound t heorem, on the other ha nd, using the theory by IChen fc Han 
1 19881) and by [Chen & Liu (119901) . Holsa pple (120081) developed a practical technique. His 
theorem guarantees that on the assumption that the velocity field is linear, for a given volume, 
the yield condition of the averaged stress over it is equal to an upper bound condition for 
the actual failure. The averaged principal stress over a given volume 14 is given as 


(Tjk 


1 

V k 



(18) 
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Here, we average the stresses over the surface shell and the internal core. Using Equation 
(USD, we calculate the averaged stress components of the surface shell as 


°sl 


L R ' 
a S 2 = -«1 — 


Os3 — — 


Rl 

Rl — R\ _|_ ^ i 


K. 1 


Rb 


5 1 

Rl 


RV 


Rl 


5 1 


RV 


and those of the internal core as 


°cl 


= <Tc 2 = h (1 


o\ 

fl") + yu 




(19) 

( 20 ) 

( 21 ) 

( 22 ) 


The details of these equations are in Appendix |A] To verify these stress comp onents, we 
show that they recover the averaged stresses over the whole volume obtained by Ho ls apple 
(120071 1. If Rb = 0, the averaged stress over the surf ace shell times np 2 GR 2 is identical to 
Equation (4.12) with a — b — c — R in Holsapple ( 2007 1. If Rb = 1, Equation (4.12) is 
recovered again since this case is reduced to aggregates with a homogeneous interior. This 
consistency allows us to compare the upper bound conditions of the surface shell and the 
internal core. The yield conditions of these stress components provide the lowest estimated 
value of the actual cohesion, later known as the lowest cohesion. Similar to Equation (lT7lh 
this is described as 


Yk = -(ofci - cr k3 ) sec (p + + vm) tan0. 


(23) 


In Section 14.11 we discuss the upper and lower bound conditions for the surface shell 
and the upper bound condition for the internal core. For the internal core, since the surface 
shedding mode is described by a core that does not reach the yield, the lower bound condition 
for the internal core does not provide constraints on it. Therefore, we do not consider this 
condition in the present problem. 


3. Numerical Modeling of Surface Granular Flow 


The simulation program that is used for this research applies a SSDEM to simulate a 
self-g ravitating granular aggregate (I Cunda 1 ll 1 1971 1: ICundall fe Hartlll992l; ISanchez fe Scheeres 
20 111). The particles, modeled as spheres that follow a predetermined size distribution, 
interact through a soft-repulsive potential when in contact. This method considers that 
two particles are in contact w hen t he y over l ap. Wh en t his happens, normal and tangential 
contact forces are calculated (IHerrmann fe Ludind l l998l ). 
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The calculation of the normal forces between colliding particles is modeled by a hertzian 
spring and a dashpot. The elastic force is modeled as 


f e = k n ^ 3/2 h, 

the damping force as 

td = -7n^n, 

and the cohesive force between the particles is calculated as 


t = -2tt 


r 2 2 

' 1 ' 2 - 
9 , 9 ^yy^-12 

r{ + r2 


(24) 


(25) 


(26) 


where rq and r 2 are the radii of the two particles in contact, a yy is the tensile strength of 
this cont act, w hich is given b y a co hesive matrix formed by the (non simulated) interstitial 
regolith flSanchez fe Scheeresl 2014), and fq 2 is the branch vector between the centers of 
these two particles. Then, the total normal force is calculated as f n = f e + f c + f^. In these 
equations, k n is the elastic constant, £ is the overlap of the particles, q n is the damping 
constant (related to the dashpot), £ is the rate of deformation and n is the normal vector 
on the surfaces of the colliding particles. As the simulation only deals with spheres, iq 2 has 
the direction as n. This dashpot models the energy dissipation that occurs during a real 
collision. 

The tangential component of the contact force models surface friction statically and 
dynamically. This is calculated by placing a line ar sp rin g attached to bot h parti cl es at 


the contact point at the beginning of the collision flHerrmann fe Luding 19981 ; iSi lbert et al. 


20011 ) and by producing a restoring frictional force f f . The magnitude of the elongation 
of this tangential spring is truncated in order to satisfy the local Coulomb yield criterion 
|ft| < 


Rolling friction (jAi et al . 2 0111) has also been implemented in order to mimic the be¬ 
havior of aggregates formed by non-spherical grains. That is, particles are subjected to a 
torque that opposes the relative rotation of any two particles in contact. This torque, simi¬ 
lar to surface-surface friction, is implemented as linearly dependent on the relative angular 
displacement of any two particles in contact and has a limiting value of: 


M™ = n r R r |f n | 


(27) 


where /i r is the coefficient of rolling resistance, R r = rir 2 /(ri + r 2 ) is the rolling radius 
and rq and r 2 are the radii of the two particles in contact. This allows our simulations 
to reach fricti on ang l es o f up to ~ 35° as evaluated by the Druker-Prager yield criterion 
flSanchez fc Scheeresl 2012). This value is typical of geological aggregates, though friction 


















angles of ~40° are not rare. This implementation of rolling friction is similar to that of 
surface friction, but is instead related to the relative angular displacement. 

To obtain the initial configuration of the aggregates, we leave the particles (initially 
cohesionless and frictionless) to coalesce only under the influence of their mutual gravitational 
interactions. Then, the particles are encapsulated inside a perfectly solid sphere that shrinks 
to a predetermined size. This size has been pre-determined from previous simulations so that 
the particles are forced into a close-to-perfect spherical shape. When the set-up procedure is 
finished, friction and cohesive forces are applied and the simulation is ready to start. Once 
the simulation starts, the aggregate is spun up in small, discrete increments of a normalized 
spin rate of 5.5 x 10~ 3 every 3000 s around the spin axis and the disruption process is 
observed. This time step showed to be long enough to allow the aggregates to safely reshape 
before the next spin-up event. The aggregates rotate in the -z direction to facilitate the 
observation of the reshaping and disruption process through the graphical interface of the 
simulation code. 

For this research, to observe the effect of an internal core in a self-gravitating aggregate, 
we have chosen to have cores whose radii are 0.5, 0.6, 0.7, 0.8 and 0.9. The particles forming 
the core are subjected to cohesive forces that are five times as strong as those felt by the 
particles of the shell. All the aggregates used in the simulations are formed by 3000 perfectly 
spherical particles with normalized radii between 0.025 - 0.035. The normalized parameters 
of the model are as follow: k n = 6.1 x 10', k t = 2.1 x 10 6 , y n = 67, /i = 0.5 and /i r = 0.8. 


4. Analysis of Surface Shedding 
4.1. Failure Condition of the Internal Structure 


With the analytical model developed above, we now investigate the effect of the internal 
core on the surface failure and the surface shedding condition. The fr iction ang l e is fixed 
at 35 °, a mean value of typical friction angles of geological materials ( Lambe fe W h itman 


1969t) . For a choice of Poisson’s ratio, we account for the earlier studies by iHolsapplel (120081 1 


and bv iHirabavashi fe Scheeres ( 20141 1. They confirmed that Poisson’s ratio do not critically 
change the failure conditiorQ. Here, to take into account compressibility of a soil material, 
we will define Poisson’s ratio as 0.25. With these fixed values, we derive the upper and lower 
bound conditions for the surface shell and the upper bound condition for the internal core. 
It is noted that since the surface shedding mode is characterized by an internal core being 


1 More specifically, they compared two Poisson’s ratios, 0.2 and 0.33, and confirmed that the critical spin 
conditions only changed within 2%. 
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below the yield and a failing surface shell, the lower bound condition for the internal core 
does not give any constraint on cohesion. 


Figure [2] shows the upper and lower conditions for the normalized cohesion as a function 
of the normalized core radius. The spin rate is constant on each line. We have also plotted 
the 1.15 spin rate, at which a particle resting on the surface can be gently lofted, and the 
1.63 spin rate, at which a particle escapes from the gravity of the test body. The former and 
latter spin rates are given by the dotted and dashed lines, respectively. Figure 2(a) gives the 
failure condition of the surface shell. The bold lines define the lowest cohesion at a given 
spin rate, while the narrow lines provide the highest cohesion at that spin rate. The actual 
failure spin rate should be between these two lines. At a constant spin rate, as the core 
radius increases, the highest and lowest cohesion decrease. On the other hand, though the 
surface shell should fail structurally, the internal core should be below the yield (see Figure 


2(b)). The bold lines describe the lowest cohesion of the internal core at a given spin rate, 


while the narrower lines shows the highest cohesion of the surface shell at that spin rate. 
These results show that at a constant spin rate, the lowest cohesion that the internal core 
needs to remain stable is higher than the cohesion that the surface shell requires to be able 
to fail. This means that if the cohesion of the internal core is not high enough, it will tend 
to fail before the surface. In a homogeneous body, failure should start therefore at its center 
and propagate towards the surface. This scen ario corresponds to the possible failure mode 
of asteroid 1950 DA flHirabayashi fc Scheeres! 120 151). 


Figures 2(a) and 2(b) imply that if surface shedding is observed, details of the process 
could provide information about their internal structure. Given a constant spin period, we 
could give constraints on cohesive strength. Since the actual cohesion should be enclosed 
by the lowest and highest cohesion curves, it should be possible to obtain the range of 
cohesion based on a possible range of the core radius. In general, the core radius is not 
well determined; however, it is still possible to give a constraint on it. For example, let us 
consider a test body that experiences surface shedding at a spin rate of 1.3. For the surface 
shell (Figure 2(a)), taking the minimum value of the lowest cohesion (0.04 at a core radius 
of 1.0) and the maximum value of the highest cohesion (0.3 at that of 0), we roughly obtain 
the range of cohesion of the surface shell as 0.04 - 0.3. On the other hand, for the internal 
core (Figure 2(b)), we roughly give that of the internal core as > 0.14. 


4.2. Dynamics of Failure Mode 

This section investigates dynamics of failure modes by changing the size of the inter¬ 
nal core. The initial conditions of the friction angle and the porosity are given based on 
experimental and observational evidences. The friction angle considered in the numerical 
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simulation is 35°, which is consistent with the analytical model. The porosity i s calc ul ated 


to be 34%, similar to that of the lunar regolith observed by Apollo 15 mission flMitchell et al. 


1972 1). This is calculated by dividing the total mass of the particles into the volume obtained 


from semi-axes of a dynamically-equivalent equal-volume ellipsoid (DEEVE) with the same 
moments of inertia as the aggregate. 


The normalized cohesion of the internal core and the surface shell is 0.5 and 0.1, respec- 
tivelyB To obtain these values, we conduct the following numerical settings; given the tensile 
strength between the particles of the shell, <j t s, the tensile strength between the particles in 
the core is described as a tc = 5 a ts and that between the particles in the interface is obtained 
as a t i — 1.5ot s . These values of cohesion were chosen as previous simulations of homogeneous 
aggregates showed their critical spin rates and disruption modes to be very different; a 0.1 
cohesion homogeneous aggregate fails at 1.13 (Figure [5]) and a 0.5 cohesion aggregate breaks 
at 1.44 (Figured]). These figures are explained later. Also, as seen in Figure [2] from the 
analytical model, known are the spin rates at which homogeneous aggregates would deform 
and disrupt, and so the ratio of cohesive forces warrants that the core would not fail before 
the surface. 


It is important to note that the yield condition and irreversible flow in the SSDEM are 
consistent w ith the_Drucker-Prager yield criterion and an associated flow rule, respectively 
( Sanchez fe Scheeres 201 2). while the analytical mode considers the Mohr-Coulomb yield 
criterion and an associated flow rule. In the present study, we use these different yield 
criteria to show how they are consistent with each other in the present problem. 


We first compare the failure conditions obtained from the SSDEM and the analytical 
models. Figure [3] shows the spin-up profile for the cases of normalized core radii of 0, 0.5, 0.6, 
0.7, 0.8 and 0.9, respectively. The peak spin rates for these cases are 1.13, 1.16, 1.18, 1.21, 
1.27 and 1.33, respectively. After experiencing their peak spin rates, the test body cannot 
reach this spin rate again. This comes from the fact that the shape has already changed 
irreversibly, causing the test body to be more susceptible to failure. An increase in spin rate 
induces further deformation with the consequent slow-down in spin rate to conserve angular 
momentum. Figure 0] shows the spin rate at a constant surface cohesion of 0.10 as a function 
of the core radius. The solid lines are obtained by the analytical model, while the dotted 
line is given by the SSDEM. The simulation results are located between the lower and upper 
bound conditions calculated by the analytical model, showing the consistency of these two 
approaches 


Next, we observe the deformation of the test body that occurs right after the rapid 


2 Since cohesion is given in units of pressure, i.e., stress components, we normalize it by tt p 2 GK 2 . 
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spin-rate drop to investigate the initial process of its failure mode. Figure [5] shows the 
deformation mode of the test body with a normalized core radius of 0. Figure 5(a) shows 
the 3-dimensional depiction of the aggregate after failure. The color shows the latitude 
of the test body. Figures 5(b) through 5(d) show the deformation vectors over the thin 
slices that include the center of mass and which are perpendicular to the x through z axes, 
respectively. Figures 5(b) and 5(c) (views from +x and +y) show how the poles of the 
aggregate push inwards towards the center, helping in the formation of an equatorial bulge. 
Figure 5(d) (view from +z) shows how the aggregate deforms ou twards from th e center. 
This deformation mode corresponds to the deformation of 1950 DA flHirabavashi fe Scheeres 


20151 ) and could lead to its catastrophic disruption of the internal core. A possible disruption 


mode for this case will be discussed in Section [5l 


As seen in Figures El and El on the other hand, the effect of the internal core on the 
deformation is obvious. The deformation of the internal core is minimal compared to the 
homogeneous (zero-core) case, while the surface shell fails structural^. For the case of a 
normalized radius of 0.5 (Figure El), after the surface shell fails, the particles there flow over 
the surface to reach the equatorial region. It is found that the failure region is not global, 
but local. For the case of a normalized radius of 0.9 (Figure El), the failure region becomes 
smaller than that for the case of a normalized radius of 0.5 and mainly occurs near the 
equatorial region. 


Now we analyze Figures El through El more closely to understand the effect of the internal 
core. These cases show the following three characteristics. First, the flow of particles is 
asymmetric with respect to the rotation axis and is symmetric with respect to the the 
equatorial plane, causing a ridge-like shape. Second, the region with the greatest deformation 
is flat and is diametrically opposed to the region that is least deformed and that is still very 
round. Third, points on the equatorial region between these two extremes show intermediate 
degrees of deformation/stretching (oblateness). These characteristics become less and less 
pronounced as the core grows in size. For a radius of 0.9, there is no deformation towards 
an oblate shape and the particles of the equator are simply ejected from the main body. 
The deformation of the test body happens only due to the failure of the surface shell, so the 
thicker it is, the greater the deformation. This implies that for_the case of a normalized core 
radius of 0, the deformation is global (jSanchez & Scheeres 2012). 


Our simulations, in no case show a deformation towards a symmetric oblate body and 
instead favor a wedge-like shape. As homogeneous as we would like our aggregates to be, by 


3 The maximum deformation of the core occurs in the equatorial plane. For a normalized core of 0.5, the 
deformation is up to 3.5% in size, which is negligible compared to the surface deformation. As the core size 
becomes larger, the deformation becomes smaller. 
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construction they will present internal weak points where the body (as any real aggregate) can 
fail more easily. This is the reason for local nature of the region of greater deformation. Once 
the aggregate deforms in one region, its internal stress is diminished as the angular velocity 
has to decrease to conserve angular momentum. In this paper we have only focused on the 
initial phase on deformation so that the dynamics of surface shedding can be understood. 
However, future efforts will be directed towards understanding of its stable shape. 

The simulations also show the dynamics of surface shedding. In Figure U the lower 
dashed line represents the condition at which the particles are just lofted from the surface, 
i.e., 1.15. If the spin rate is faster than this line and the surface shell already fails, surface 
shedding should occur. The upper dashed line represents the condition at which particles 
escape to infinity, i.e., 1.63. According to this figure, if the core radius is 0, particles cannot 
be shed. For the cases of core radii of 0.5 and larger, particles are lofted. However, since 
all the cases are below the upper dashed line, particles cannot reach escape velocity. They 
would either stay in orbit around the primary body or come back to its surface. Potentially, 
an absolute increase in the cohesion of the shell and interior would enable shed particles to 
immediately escape upon their being shed. 


Discussion 


In this paper, we have proposed a possible structure for a rubble pile asteroid that 
experiences surface shedding. It_is im portan t to consider what we know about failure modes 
from ground observations. Holsapple (120101) found that if the internal structure is uniformly 
distributed, the original body completely deforms before surface shedding. This predicts 
that a homogeneous spheroid would never have surface shedding, b ut co uld experience a 
catastrophic breakup at high enough cohesion. 1 Sanchez & Schee resl (l2014al Jbll simulated the 
zero core case for 560000s, which is 10000s longer than the initial failure time (~ 460000s). 
The cohesion was fixed at 0.50, which is the same cohesion used for the internal core in the 
previous sections. Figure [8] shows a snapshot of the failure mode at 560000s. The test body 
catastrophically breaks up into multi-components. This mode is clearly distinguished from 
surface shedding. The surface shedding mode originates from the failure of the surface shell, 
causing granular flow. 

Several asteroids recently observed are now considered to have undergone surface shed¬ 
ding. A stero i d P/2013 P5 (PA NSTARRS) was reported to experience its episodic dust 
ejection ( Jewitt et al. 2013 . 2015h . This asteroid is orbiting near the inner edge of the aster¬ 
oid belt in the vicinity of the Flora family of S' -type asteroids. It implies that sublimation of 
water ice is unlikely to occur on this asteroid (Hewitt et al.ll20131 ). The size of dust particles 
was up to a few millimeters, negligibly smaller than the nucleus, and the ejection velocity 
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was less than 1 m/s (IJewitt et al.ll2015h . From these observational results, the episodic event 
of this asteroid is currently considered to come from surface shedding. Asteroids 133P/Elst- 
Pizarro and (62412) 2000 SY178 have relatively rapid spin periods of 3.47 hr and 3.33 hr, 
respectively, and had experienced mass ejection. Although the mass ejection might have 
been driven by the seasonal activity, i. e., sublimat i on of water ice, rotational instability i s 
also considered to help their activity flHsieh et al.l 12004 l2010l; ISheppard fe Truiillol 120151) . 
Some remarkable observations related to their activity include: (1) the size of ejecta is small, 
(2) the mass ejection is negligible compared to the nucleus and (3) the ejection velocity is 
relatively small. These facts contradict the catastrophic breakup of the coreless spheroid 
observed in Figure [HJ but is consistent with our model of an internal core. It implies that 
these active asteroids may have internal cores which prevent them from breaking up catas¬ 
trophically. 

Additionally, simulations allowed us to observe that the failure due to rotational defor¬ 
mation has a local nature. The local deformation of the surface shell is dependent on its 
structure, especially on the spatial distribution of particles which is translated in the spatial 
variation of porosity and the possible presence of fractures. Even though the settling pro¬ 
cedure, devoided of cohesive and frictional forces, encourages structural homogeneity, this is 
not perfect and, at high enough spin rates, granular flows always occur locally. We consider 
this to be the initial deformation mode. Though long-term deformation modes are not the 
focus of this paper, there are two possible scenarios for them. In the first scenario, material 
flows occur globally, causing the formation of a unifor m equ atorial ridge. This scenario cor¬ 
responds to the result by obtained Walsh et ah (120081 . 2012). In the second scenario, a local 
deformation leads to further structural bias, causing the further emergence of local material 
flows. We will conduct further investigations for this problem in the future. 

The formation process of a strong core is beyond our scope; thus, we do not quanti¬ 
tatively analyze it in this paper. However, there may be two likely scenarios. The first 
scenario is that at an accretion process after a catastrophic impact, due to the difference of 
energy dissipation, large boulders and rocks could accumulate first and then small particles 
might cover these large objects. Such a formation process could produce bodies with layers 
of different strength similar to the model we have here analyzed (Durda and Walsh, 2015 
personal communication). Note that since m echan ical strength depends on many different 
parameters such as bulk density flHirabavashil 20141 ). friction angle and cohesion, we cannot 
conclude that such a structure gives a strong core. We leave this problem for future works. 
The second scenario is based on the assumption that the distribution of loose materials 
is uniform initially. This condition leads to compression of loose particles in the interior. 
When such particles irreversibly mo ve out of their original locations, the volume tric strain 
decreases as compression increases flLambe fe Whitman! 19691 : Nedderman 1992 ). Because 
of the self-gravity, the interior experiences larger compression than the surface. As a result, 
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parti cles in the interi or ar e highly packed. Since higher compaction gives higher cohesion 
flLambe fe Whitman! 1 9691 ). this scenario could produce a strong core in an asteroid. 


F inally , we compare the HSDEM bv I Walsh et al.l (120081.120121) with our model. IWalsh et al 


(l2008l . 2012) observed surface shedding from a rubble pile body with a homogeneous hexago¬ 
nal closed packing (HCP) structure^, while we concluded that a body should have an internal 
core to have sur face s hedding. This different result comes from the fact that the simulations 
carried out by Walsh et al. ( 20081 . 120 12!) and our model may have different deformation fea¬ 
tures. If particles are arranged to form a crystalline lattice, each particle disl ocation is 


highly controlled by its closest neighbors and the deformation is not isotropic (lAshbv et ah 


20141 ). However, even for such a configuration, particles sitting on the surface have open 
spaces where they can move freely; thus, their motion is less restricted than if they were in 
the interior. Consequently, the initial HCP configuration could cause anisotropic flows and 
provide t heir aggre gates with a strong core and a weak surface shell implicitly. Note that 
Walsh et al. (2012) also investigated a near-fluid case that was modeled by the bi-modal 
size distribution (a randomly packed aggregate) and obtained a pancake shape, which might 
indicate an isotropic deformation mode. 


6. Conclusion 

This paper quantitatively analyzed the effect of the internal core of a rubble pile spheroid 
on the mechanism of surface shedding. On the assumption that the friction angle and the 
bulk density are uniform over the whole volume, we introduced a simple surface shedding 
model. It is composed of two different cohesion layers: the internal core and the surface 
shell. The analytical method was constructed by the lower and upper bound theorems 
in limit analysis, while the numerical method was based on the SSDEM by Sanchez and 
Scheeres. The failure conditions derived by the analytical model were consistent with those 
by the SSDEM. The results showed that an aggregate with homogeneous cohesion would 
always fails internally before it fails superficially. To have surface shedding, therefore, the 
progenitor body needed a strong core so that only the surface region failed structurally. The 
SSDEM also exposed the dynamics of the initial failure phase. The primary result was that 
for all the cases, a rubble pile spheroid would fail locally through a severe material flow 
on its surface. This came from the fact that since the initial body configuration was not 
perfectly uniform, which we believe is more realistic, some surface regions are susceptible to 
fail earlier than others. Further research would be able to provide a better understanding of 
the formation of unique shapes of asteroids and the mechanisms that originate the activity 


Whey defined several initial packing configurations to control the friction angle. 
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observed in active asteroids. 


The authors wish to thank Dr. Kevin Walsh and Dr. Dan Durda at SwRI Boulder and 
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A. Averaged Stresses over the Core and the Surface Shell 


The averaged stresses over the internal core and the surface shell are given by the inte¬ 
grations of Equations (J2J) through (|7J) over these volumes. Since the shear stress components 
given by Equations (J5J) through (|7J) are trigonometric functions, they become zero by the 
integration. On the other hand, the normal stress components, Equations (J2J) through (jTJ), 
have non-zero values. Using the indices notations defined above, we obtain the averaged 
normal stresses for the surface shell as 
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and those for the internal core as 
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Note that an = <T 2 S and o\ c = a^c- Since the shear stress components are zero, these 
components are considered to be the principal stress components. 
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Fig. 1.— Sphere model with the internal core. The sphere is assumed to be spinning 
constantly along the z axis. The normalized radii of the sphere and the internal core are 
given as 1 and Ri,, respectively. 
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(a) 



Normalized Core Radius 


(b) 


Fig. 2.— Normalized cohesion as a function of the normalized core radius. The numbers are 
the normalized spin states, which are constant on each line. Figure 2(a) shows the highest 
(bold lines) and lowest (narrow lines) cohesion of the surface shell. Figure 2(b) describes 
the highest cohesion of the internal core (bold lines) and the lowest cohesion of the surface 
shell (narrow lines). The dotted lines indicate the spin rate at which a particle is just lofted, 
~ 1.15, while the dashed lines describe that at which a particle escapes to infinity, ~ 1.63. 
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Fig. 3.— Time evolution of the angular velocity of the test body. The colors indicate 
different core radii; the red, the green, blue, magenta, light blue and yellow show normalized 
core radii of 0, 0.5, 0.6, 0.7, 0.8 and 0.9, respectively. The angular velocity of the test body 
is incremented in steps of a normalized spin rate of 5.5 x 10 -3 every 3000 s. 
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Fig. 4.— Comparison between the analytical model (the solid lines) and the SSDEM (the 
dotted line). On these curves, the cohesion is hxed at 0.10. The SSDEM result is perfectly 
enclosed by the upper and lower bound conditions derived by the analytical model. The 
lower dashed line indicates the condition at which particles are just shed from the surface, 
while the upper dashed line describes that when particles escape. 
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Fig. 5.— Deformation of the test body with a normalized core radius of 0. Thus, the 
normalized cohesion of the whole body is 0.1. The normalized critical spin rate is 1.13. 
Figure 5(a) shows the 3-dimensional view. The color shows the latitude of the test body. 
Figures 5(b) through 5(d) show the displacement vectors of the particles contained in thin 
slices perpendicular to the x, y and z axes as seen from +x, +y and +z respectively. In these 
pictures, the blue dots are the original locations of the particles. 
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Fig. 6.— Deformation of the test body with a normalized core radius of 0.5. The normalized 
cohesion of the surface shell and that of the internal core are 0.1 and 0.5, respectively. The 
normalized critical spin rate is 1.16. The definitions of the figures are the same as Figure [5] 
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Fig. 7.— Deformation of the test body with a normalized core radius of 0.9. The normalized 
cohesion of the surface shell and that of the internal core are 0.1 and 0.5, respectively. The 
normalized critical spin rate is 1.33. The definitions of the figures are the same as Figure [5l 
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t=560000?0004^-— 

0Mi=-0.000030rad/s 
OMf=-0.001296rad/s 


Fig. 8.— Example of the final failure mode of a zero-core spheroid (jSanchez fc Scheeres 
2014a|jbl l. The normalized cohesion of the whole body is 0.5, and the normalized critical spin 


rate is 1.44. The simulation time is 560000s, which is 10000s longer than the time when 
the initial failure mode is observed (~460000s). The progenitor body is completely breaking 
into several large components. 
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